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We have measured the form factor ratios ry ~ V{0)/ Ai{0) and r2 — ^2(0)/y4i(0) for the decay 
^ (j) vt, (f) ^ K^ K~ , using data from charm hadroproduction experiment E791 at Fermilab. 
pH I Results are based on 144 signal and 22 background events in the electron channel and 127 signal 

and 34 background events in the muon channel. We combine the measurements from both lepton 
channels to obtain ry = 2.27 ± 0.35 ± 0.22 and r2 = 1.57 ± 0.25 ± 0.19. 
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Semileptonic decays of mesons containing one heavy 
quark are described in simplest terms by a spectator 
model in which the heavy quark decays semilcptonically 
while the light spectator antiquark passes from the par- 
ent meson to the final-state meson. These decays are still 
relatively simple when one considers hadrons rather than 
quarks; the leptonic and hadronic currents factorize, and 
all strong-interaction effects can be parametrized by a 
few form factors 

If SU(3) flavor symmetry is approximately valid, re- 
placing a spectator d quark by a spectator s quark 
should have little effect on the form factors. We have in- 
vestigated this spectator replacement experimentally in 
charm meson decay, and report here new measurements 
by Fermilab experiment E791 of form factor ratios for 
D+ (jje+Ve and D+ (jjfJ-'^i^f,, with K+K^. 
These Ds results are then combined and compared to 
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recent high-statistics form-factor ratio results for 

D"*" K i'^ i/g measured in the same experiment, and 
also to theoretical predictions. (Charge-conjugate states 
are also included by inference throughout this letter.) 

In addition to testing SU(3) symmetry in semileptonic 
decays, ratios of form factors in D+ ^ (j)£^ I'g are needed 
for tests of factorization in Ds hadronic decays [Q . 

With a vector meson in the final state, there are four 
form factors, V{q'^), Ai{q'^), A2(g^) and ^3(9^), which 
are functions of the Lorentz-invariant momentum trans- 
fer squared Q| . The limited size of current data samples 
precludes precise measurement of the g^-dependence of 
form factors, and thus we assume a nearest-pole dom- 
inance model: F{q^) — F(0)/(1 — <f/^'^oie) where 
rupoie = mv ~ 2.1 GeV/c^ for the vector form fac- 
tor V, and TTipoie = raA = 2.5 GeV/c^ for the axial- 
vector form factors jsj. Because Ai{q^) appears among 
the coefficients of every term in the differential decay 
rate, it is customary to factor out Ai[Q) and to mea- 
sure the ratios ry = V{Q)/Ai{Q), ra = A2(0)/Ai(0) and 
r3 = A^{Q)/Ai{0) which are independent of the absolute 
normalization of the data. 

The form factor A3, which probes the spin-0 compo- 
nent of the virtual W , is unobservable in the limit of 
vanishing lepton mass, and has been found to be small 

in £)+ K /i+i^p [||. However, the q^ dependence of 
the muon channel is more complex than that of the elec- 
tron channel even if Aj, vanishes. The decay rate at low 
q^ is suppressed by an overall factor (1 — m'jjq'^)'^, and 
additional spin-flip amplitudes, suppressed by an over- 
all factor of rn^^/q^ relative to the spin non-flip ampli- 
tudes, contribute to the differential decay rate at small 
q^ . E791 is the first experiment to measure form factor 
ratios for semileptonic decays in both electron and 
muon modes. Obtaining consistent results for the two 
lepton channels gives confidence that these q^ differences 
are adequately understood. 

Fermilab E791 is a fixed-target charm hadroproduc- 
tion experiment Charm particles were produced in 
collisions of a 500GeV/c beam with five thin targets, 
one platinum and four diamond. About 2 x IQ^'^ events 
were recorded during the 1991-1992 Fermilab fixed-target 
run. The tracking system consisted of 23 planes of silicon 
microstrip detectors, 45 planes of drift and proportional 
wire chambers, and two large-aperture dipole magnets. 
Hadron identification is based on the information from 
two multicell Cerenkov counters that provide good dis- 
crimination between kaons and pions in the momentum 
range 6 — 36 GeV/c. In this momentum range, the prob- 
ability of misidentifying a pion as a kaon depends on 
momentum but does not exceed 5%. 

Electron candidates are identified by electromagnetic 
shower shape, the match between calorimeter energy and 
tracking momentum, and agreement between calorime- 
ter and tracking position measurements. Muon can- 
didates are identified by a single plane of scintillator 
strips, oriented horizontally, located behind an equiva- 
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lent of 2.4 meters of iron (comprised of the calorime- 
ters and one meter of bulk steel shielding). The an- 
gular acceptance of the scintillator plane was about 
±62 mrad x ±48 mrad (horizontally and vertically, re- 
spectively), which was somewhat smaller than that of 
the rest of the spectrometer for tracks that went through 
both magnets (« ±100 mrad x ±64 mrad). The verti- 
cal position of a scintillator hit is determined from the 
strip's vertical position, and the horizontal position of a 
hit from timing information. 

The event selection criteria used for this analysis are 
the same for both electron and muon channels except 
for those related to lepton identification. Events arc se- 
lected if they contain an acceptable decay vertex deter- 
mined by the intersection point of three tracks that have 
been identified as a lepton and two kaons. The longitudi- 
nal separation between this candidate decay vertex and 
the reconstructed production vertex is required to be at 
least 11 times the estimated error on the separation, and 
the candidate decay vertex must occur outside the target 
material. 

We require electron candidates to have a minimum mo- 
mentum of 5 GeV/c, which gives an electron identifica- 
tion probability of about 70% and a probabability for a 
pion to be misidentified as an electron of 1-2%. We re- 
tain only muon candidates with momenta greater than 8 
GeV/c to reduce the contamination from hadron decays 
in fiight. With this momentum restriction, the muon 
identification efficiency is about 85%, and the probabil- 
ity for a hadron to be misidentified as a muon is about 
3%. 

Events in which the lepton candidate has the same 
charge as the candidate decay are assigned to a "right- 
sign" (RS) sample. This sample contains both signal 
events and a small amount of background arising from 
reconstruction errors, including misreconstructcd charm, 
and possible kinematic reflections from other charm de- 
cay channels. Events in which the lepton candidate and 
the candidate decay have opposite charge are assigned 
to a "wrong-sign" (WS) channel, which contains purely 
background. Because WS events have two hadrons of the 
same charge, this sample is dominated by reconstriiction 
background. We expect the amount of such background 
in the RS sample to be about twice the number of events 
in the WS sample, since two tracks rather than one can 
be misidentified as the candidate lepton. 

To cxchide fecdthrough from the hadronic decays 
D+ K-K+TT+ and D+ -> K-K+Tr+, we exclude 
events in which the invariant mass of the three charged 
particles (with the lepton candidate interpreted as a 
pion) is consistent with the £>+ or Df mass. To exclude 
feedthrough from the more copious semileptonic decay 
D+ ^ K e+Pe, we exclude events in which the invari- 
ant mass of the two hadrons, with K~ and tt"*" masses 

^ Q 

assigned, is consistent with K (892). 

Two kinematic signatures display the D'^ (l)£'^i'e 
signal: a Dg peak in the minimum parent mass M^i^i 



3 



and a 4> peak in K~ invariant mass. Mmm is de- 
fined as the invariant mass of KKtui wlien the neutrino 
momentum component along the D'^ direction of flight 
is ignored. The Mmin distributions for the electron and 
muon channels are shown in the top left plots of Figs. 
1 and 2, respectively; Dg peaks are evident in the RS 
events. We retain events with Mmin in the range 1.7-2.1 
GeV/c^, as indicated by arrows in the figures. The dis- 
tribution of KK invariant mass for the retained events 
is shown in the top right plots of Figs. 1 and 2; the RS 
sample exhibits a clear </> mass peak. Candidates with 
1.014 < Mkk < 1-026 GeV/c^ were retained, yielding 
final data samples of 166 right-sign and 11 wrong-sign 
events for the electron channel and 161 right-sign and 17 
wrong-sign events for the muon channel. 

A Monte Carlo calculation simulated the response of 
the E791 detector to semileptonic decays. A sam- 
ple of events was generated according to the differential 
decay rate (Eq. 22 in Ref. |^), with form factor ratios 
r2 — 0.78, ry = 1.74, and — 0. The same selection 
criteria were applied to the Monte Carlo events as to 
real data. Out of 4.5 million events generated in the two 
modes, 13640 decays passed all cuts. The bottom plots 
in Figs. 1 and 2 compare the background-subtracted 
(RS-2xWS) (j) signals to the Monte Carlo K+R- spec- 
tra from — !• vi decays. The agreement between 
the two distributions indicates that the background un- 
der the (f) peak is adequately described by twice the WS 
events. 

The differential decay rate |0] is expressed in terms 
of four independent kinematic variables: the square of 
the momentum transfer (q^), the polar angle Oy between 
the and in the rest frame, the polar angle Oi 
between the vi and D'^ in the rest frame, and the 
azimuthal angle % between the (p and W'^ decay planes 
in the rest frame. Our definition of the polar angle Og 
is related to the definition used in Ref. ^ by 6*^ ^ n — 9e. 
We do not use x in this analysis because our data sample 
is not large enough to profit from the small increase in 
sensitivity that x can provide. 

Semileptonic decays cannot be fully reconstructed due 
to the undetected neutrino. With the available in- 
formation about the D+ direction of flight and the 
charged daughter-particle momenta, the neutrino mo- 
mentum (and all the decay's kinematic variables) can be 
determined up to a two-fold ambiguity if the parent mass 
is constrained. Monte Carlo studies show that the differ- 
ential decay rate is more accurately determined, on av- 
erage, if it is calculated with the solution corresponding 
to the lower laboratory-frame neutrino momentum. 

To extract the form factor ratios, the distribution of 
data points in the three-dimensional kinematic variable 
space is fit to the full expression for the differential de- 
cay rate including lepton mass terms, but with A3 set 
to zero. We use an unbinned maximum-likelihood fitting 
technique in which the likelihood function is computed 
from the density of Monte Carlo events near each data 
event in the three-dimensional space of kinematic vari- 
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ables; the same Monte Carlo events are reweighted to 
provide different trial values of form factor ratios. To in- 
clude background in the fit, a similar likelihood function 
based on twice the density of wrong-sign events around 
each right-sign event is used. 

With the above method, the fitted results are sub- 
ject to small systematic biases which originate from two 
sources: (a) approximate normalization of the likelihood 
function and (b) nonlinearity of the decay rate within the 
volume centered on the data point. These systematic bi- 
ases of the fitted parameters were determined from Monte 
Carlo studies: Sry = -HO. 18 ±0.10 and (5r2 = -0.12±0.05 
for the electron channel, and dry = -1-0.14 ± 0.09 and 
Sr2 = —0.10 ± 0.06 for the muon channel. After correc- 
tion for these biases by subtracting these (5r's from the 
measured values, the final form factor ratios and their 
statistical errors are ry = 2.24±0.47 and r2 = 1.64±0.34 
with a correlation coefficient of —0.15 for the electron 
channel, and ry = 2.31 ± 0.54 and r2 = 1.49 ± 0.36 with 
a correlation coefhcient of —0.20 for the muon channel. 

Fig. 3 compares the background-subtracted data in the 
kinematic variables , cos Og and cos dy for the electron 
and muon channels with Monte Carlo predictions using 
our measured form factor ratios. Note that although the 
distributions for the two channels are quite different, 
due in large part to lepton-mass effects, the form factor 
ratios agree well. 

The most important sources of systematic uncertainty 
and their estimated contributions to the errors in ry and 
r2 are listed in Table |. The largest systematic effects are 
associated with small densities of Monte Carlo events in 
kinematic variable space. We reject data events (5 to- 
tal) with zero nearby Monte Carlo events in this space; 
the dependence of the form factor ratios on requiring 
more nearby Monte Carlo events was studied by vary- 
ing the minimum number from 1 to 5. Dependence of 
ry and r2 on the fineness of partitioning of the variable 
space was investigated by varying this parameter by a 
factor of three. The systematic biases in ry and r2 de- 
scribed above were determined from ensembles of Monte- 
Carlo experiments; the small uncertainties in them are 
largely statistical. Sensitivity to possible inaccuracies in 
the Monte Carlo simulation of detector response was esti- 
mated with several sets of reasonable event-selection cri- 
teria. Finally, sensitivity to uncertainty in the amount of 
background was estimated by multiplying the wrong-sign 
events by weights ranging from 0.5 to 3.5. 

The consistency within errors of our results measured 
in the electron and muon channels indicates that lep- 
ton mass effects are handled correctly in the fits. This 
agreement also supports the assumption that strong in- 
teraction effects, incorporated in the values of form factor 
ratios, do not depend on the particular leptonic de- 
cay. We therefore combine the results measured for the 
electron and muon decay modes in a weighted average: 
ry = 2.27 ± 0.35 ± 0.22 and ra = 1.57 ± 0.25 ± 0.19. The 
statistical and systematic uncertainties of the averaged 
results were determined using the general procedure de- 
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scribed in Ref. § (Eqns. 3.40 and 3.40'). The largest 
source of systematic error has a positive correlation be- 
tween the two samples. 

Table |l| and the top portion of Fig. 4 compare the 
values of the form factor ratios ry and r2 measured by 
E791 with previous experimental results The size 

of the data sample and the decay channel are listed in 
Table || for each experiment. The E791 combined re- 
sults are consistent with the previous world averages 
and have total errors 20% - 30% smaller than the world 
averages. 



Table HI and Fig. 4 (bottom) compare the combined 
E791 results for ry and r2 to published theoretical pre- 
dictions for £)+ (pi^ vi ||l^-[T5|] and to E791 measure- 
ments [^,|| for _D+ K* f+i/f. Our value of ry for 
the Ds decay is in good agreement with all quoted pre- 
dictions, and is also consistent with the E791 measure- 
ment for the decay. For r2, our error for the Ds 
decay is significantly smaller than the spread in theoreti- 
cal predictions. Furthermore, our measurement for Dg is 
more than two standard deviations higher than the E791 
_D+ result. This discrepancy in the measured value of 
r2 for Ds and _D+ decays is even larger, about 3.3 stan- 
dard deviations, when new world averages are computed 
using data from this letter and from Refs. |^ and [ p2[ : 
ra = 1.58 ± 0.25 for D+ (/)£+ v^, and ra = 0.73 ± 0.07 

for Z?+ K vg. Since the quoted theoretical papers 
predict form factor ratios equal to within 10% for the two 
decays, it is difficult for these predictions to agree with 
both Ds and D^ data. 

To summarize, we have measured the form factor ratios 
for the semileptonic decay — > (pi'^i'i in both electron 
and muon modes. For the electron mode ry = 2.24 ± 
0.47 ±0.21 and ra = 1.64 ± 0.34 ± 0.20, and for the muon 
mode ry = 2.31 ± 0.54 ± 0.26 and ra = 1.49 ± 0.36 ± 0.20. 
The combined E791 results for Z?+ ipt^iyi are ry = 
2.27±0.35±0.22andr2 = 1.57±0.25±0.19. The ry result 
is consistent with the expected SU{3) flavor symmetry 
between Ds and D^ semileptonic decays; however, the 
r2 result appears inconsistent with this expectation. 
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TABLE I. Main contributions to systematic uncertainties 
in the form-factor ratios for — > (j) e^Ve and D'^ — > (f) jj.'^Ufj,. 

Source of uncertainty ar^^ (e, fj,) ct,.2 (e, m) 

Min. # of MC points per data point 0.16,0.19 0.17,0.13 

Size of partitions in variable space 0.06,0.11 0.06,0.10 

Systematic biases in fitting technique 0.10, 0.09 0.05, 0.06 

Simulation of detector effects 0.06, 0.08 0.06, 0.07 

Number of background events 0.04, 0.05 0.03, 0.04 

Total 0.21,0.26 0.20,0.20 
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TABLE II. Comparison of E791 results for D+ (fil'^vt 
form factor ratios with previous experimental results. 



Exp. 



Events 



r-v = V(0)Mi(0) ra = A2(0)Mi(0) 



E791 
E791 
E791 



271(e + /i) 2.27 ±0.35 ± 0.22 1.57 ± 0.25 ± 0.19 
144 (e) 



2.24 ±0.47 ±0.21 
2.31 ±0.54 ±0.26 



1.64 ±0.34 ±0.20 
1.49 ±0.36 ±0.20 



CLEO 
E653 1 
E687 


Jg 

LO 
11 


1 308 (e) 
19 (m) 
90 (m) 


0.9 ±0.6 ±0.3 
2.3tj i ±0.4 
1.8 ±0.9 ±0.2 


1.4 ±0.5 ±0.3 
2.llg-^ ±0.2 
1.1 ± 0.8 ±0.1 


Prev. 


world avg. |12 


1.5 ±0.5 


1.6 ±0.4 



TABLE III. Comparison of E791 D+ ^ (p t+vt form factor 
ratios with theoretical predictions for Dt ^ (jit^ V( and with 
E791 results for D+ -» 'K*° £+ vi. 



Group 


rv 


r2 


E791 Dt (t}e+Vi 


2.27 ±0.41 


1.57 ±0.31 


BKS Pi 


2.00±0.19t;^;25 


0.78 ± 0.08t[j:i^ 


lmmsHm 


1.65 ±0.21 


0.33 ±0.33 


ISGW2 PI 


2.1 


1.3 


E791 D+ K*" i+ui 1 




1 1.87 ±0.11 


0.73 ±0.10 
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FIG. 1. Distributions of minimum parent mass Mmin and 
KK invariant mass for Df — > (pe^Ve candidate events. 
Right-sign (RS) and wrong-sign (WS) samples arc defined in 
the text. Top left: Mmin for events with KK mass in the 
range 1.014 to 1.026 GeV/c^. Top right: KK invariant mass 
for events with Mmin in the range 1.7 to 2.1 GeV/c^. Bottom: 
background-subtracted (RS— 2xWS) KK mass distribution 
(crosses) compared to Monte Carlo prediction (dashed his- 
togram) for events with Mmin in the range 1.7 to 2.1 GeV/c^ 
and KK mass in the range 1.014 to 1.026 GeV/c^ . All can- 
didates pass all the other final selection cuts. The arrows 
indicate the range of the final sample. 
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FIG. 2. Distributions of minimum parent mass Mmin and 
KK invariant mass for Df (pfX^Vn candidate events. 
Right-sign (RS) and wrong-sign (WS) samples arc defined in 
the text. Top left: Mmin for events with KK mass in the 
range 1.014 to 1.026 GeV/c^. Top right: KK invariant mass 
for events with Mmin in the range 1.7 to 2.1 GeV/c^. Bottom: 
background-subtracted (RS— 2xWS) KK mass distribution 
(crosses) compared to Monte Carlo prediction (dashed his- 
togram) for events with Mmin in the range 1.7 to 2.1 GeV/c^ 
and KK mass in the range 1.014 to 1.026 GeV/c^ . All can- 
didates pass all the other final selection cuts. The arrows 
indicate the range of the final sample. 
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FIG. 3. Comparison of kinematic variable distributions for 
background-subtracted data (crosses) with Monte Carlo pre- 
dictions (dashed histograms) for electron (left) and muon 
(right) channels. The predictions use our best-fit values for 
the form factor ratios. 
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FIG. 4. Top: Comparison of experimental measurements 
of form factor ratios rv and r2 for Df (pl'^ve- The smaller 
error bars indicate the statistical errors and the larger ones 
show the statistical and sytematic errors added in quadrature. 
Bottom: Comparison of E791 measured form factor ratios for 
Df — > (j)t^ui with theoretical predictions for Dt — » V( 
and with E791 measurements for — > K**^ vt. 
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